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Abstract 

The hadron production in e'^e~ colhsions is studied assuming that particles orig- 
inate in a hadron gas at thermal and chemical equilibrium. The parameters of the 
hadron gas are determined with a fit to the average multiplicities of various hadron 
species measured at LEP and PEP-PETRA colliders. An impressive agreement is 
found between the predictions of this model and data for almost all particles over a 
range of production rate of four orders of magnitude. The temperature values found 
at the centre of mass energies of LEP and PEP-PETRA are around 165 MeV. 
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1 Introduction 



Hadron production in e"'"e~ collisions at high energy is generally believed to be the result of a 
two-stage process: a parton shower generated by the qq pair emerging from the annihilation, 
and a fragmentation of the partons into observable hadrons. The former process is hard and 
can be described by perturbative QCD, whereas the latter is soft and not calculable with a 
perturbative approach. 

Several phenomenological models aimed at describing quantitatively the fragmentation process 
have been developed; amongst them, the most popular are those implemented in the Monte 
Carlo generators JETSET |l|] and HERWIG 0, using the concepts of string and cluster frag- 
mentation respectively. The main unsatisfactory feature of these models, as far as the relative 
production rates of the hadron species is concerned, is the large number of free parameters 
required in order to correctly reproduce experimental data. As a consequence, those models 
have a rather poor predictive power. Empirical regularity in hadron production rates have 
been observed in ref. although some non-physical assumptions were made about meson 
quantum numbers and masses. The most interesting result of those observations is that the 
very same exponential behaviour of the hadron production function is present both at LEP 
(Vs = 91.2 GeV) and PEP-PETRA = 29 ^ 35 GeV) centre of mass energies. 
We introduce a thermodynamical approach to the problem of hadronization f^, postulating 
the existence of a hadron gas in thermodynamical equilibrium before the hadrons themselves 
decouple (freezing-out) and decay giving rise to observable particles in the detector. We show 
that this model is able to fit almost all inclusive rates measured so far at LEP and PEP-PETRA 
colliders in a natural way by using only three parameters. Two of them are the basic parame- 
ters of the hadron gas description, namely its temperature T and its volume V\ the third one, 
7s is a parameter describing the partial strangeness chemical equilibrium, which has been used 
already in some analysis of hadron production in heavy- ion collisions The only required 
inputs to determine the particle yields are the mass, the spin and the quantum numbers such 
as baryon number, strangeness, charm and beauty. 

2 Hadron gas scheme 

Most hadronic events in high energy e'^e~ annihilations are two-jet events. Here we assume that 
each jet represents a hadron gas phase in complete thermodynamical equilibrium just before 
the freezing-out time and that the number of such phases in hadronic events is always two, 
that means neglecting multi-jets events. One can then describe a jet as an object defined by 
thermodynamical and mechanical quantities such as the temperature and the volume in its rest 
frame and the Lorentz-boost factor 7. As far as chemical equilibrium is concerned, we assume 
that, in general, a jet has quantum numbers related to those of the parent quark. As the two 
jets must be colourless and do not have fractional baryon numbers, sharings of one or several 
quark- ant iquark pairs should occur. It seems reasonable to assume that most pairs are of either 
uu or dd type which are the lightest ones. We allow also ss sharing in the inter-jet interaction 
and non- vanishing baryon number for each jet provided that the baryon number and strangeness 
of the whole system are zero. Morevoer, we assume that each jet keeps the charm and beauty 
of the parent quark, i.e. no heavy quark pairs are exchanged. Adequate tools to deal with 
such a problem in a statistical mechanics framework are the canonical partition functions of 
systems with internal symmetries. The partition function of a system which transforms under 
the irreducible representation of a symmetry group G can be expressed as [0, §] : 
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M{G) J ^^^^^^''^^ ^)^i9), (1) 

where ^{g) is the group measure, M{G) = J dfi{g), d„ is the dimension of and Xv{g) is the 
character of ^ namely Xui.9) = (D'^ig))- The function Z is defined as: 

Z(^) =tr(e-'^^e-'S^'QO , (2) 

where the (/)j's are the parameters of the group and Qi its generators. The exponential factor 
is the usual canonical distribution with /? as the inverse temperature and E as the energy of 
the system, both calculated in its rest frame. 

In the present case the symmetry group is f/(l)^, each U{1) corresponding to the conservation 
of baryon number N, strangeness S, charm C and beauty B respectively. The Eq. (1) becomes: 

1 /'27r/'27r/'27r/'27r 

- (2?)! i /„/„/„ '''^^''' ^m- (3) 

Q = {N, S, C, B) is a four dimensional vector with integer components representing the quan- 
tum numbers of a jet and = 02, 03, 04), each 0j being the parameter of a U{1) group. 
According to the assumptions made for the quantum numbers of the system, the partition 
function of the whole system can be written as: 

z=jzjz ^(Q)^(-Q) , (4) 

Ar=0 5=-2 

where Z{Ql) are the partition functions of the single jets and the sum runs over baryon number 
N and strangeness S. Further constraints are implicitly introduced in Eq. (4), namely a 
maximum allowed baryon number 1 and a maximum allowed strangeness 2. Moreover, in case 
of e^e~ — >■ ss events the term with 5* = is removed from the sum. These assumptions define 
what has been called correlated jet scheme |^. For a gas of Nb boson species and Np fermion 
species, the partition function of a single jet becomes: 



Z{d) = J^,J e^'^'^exp ^ log (1 - e-^^--^^^'^)-' + 



(2vr) 

Np 

+$:Elog(l + e-^^'=-''^-*)}, (5) 
j=i k 

where k labels all available states of phase space for the j^^ particle, Sk is its energy and 
q^- = {Nj, Sj, Cj, Bj). Since the total number of each hadron in the gas is not determined, all 
chemical potentials must vanish. From Eq. (4) we are able to estimate the average number per 
jet of any hadronic particle assigning each of them a fictitious fugacity Aj which multiplies the 
g-/3efc factors and deriving: 

< rii >= Xi |a,=i . (6) 

The sum over the phase space, for a continuous level density becomes: 
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where J is the spin of the particle and V the volume of the hadron gas phase. 
One is left with complicated integrals in Eq. (5), (6) but, if T ~ O(IOO) MeV, as we can argue 
from the QCD soft scale, the exponential factors are expected to be very small for all particles 
but pions, so that 

i ~ f,-^/p^+m^/T-^(li■<t> ^g^ 



are very good approximations. This corresponds to the Boltzmann limit of Fermi and Bose 
statistics. It can be shown then that the average number of particles i produced per jet is 



^ ' Z(Q) ' 



whereas for pions 



being 



< rii >= ,^ ,„ I d p — , , 



(10) 



z^ = (2J. + 1)^ / d'p e-V5^/^ = (2J, + l)^,rn^K,C^) ■ (H) 
K2 is the modified Bessel's function of order 2. 

The factor Z(Q — qi)/Z(Q) suppresses or enhances the thermal production rate of the particle 
according to its quantum numbers and the quantum numbers of the jet (see also Fig. 1). 



3 Fit of LEP and PEP-PETRA data 

The volume and the temperature of the two phases at the decoupling time are determined by 
a fit to the available data on hadron inclusive production. The volume and the temperature 
are assumed to be equal for both jets and also equal for different kinds of primary quarks. A 
parameter 7^ < 1 is introduced in order to take into account a non-complete strange chemical 
equilibrium; if a hadron contains n strange valence quarks, its production rate is reduced by a 
factor 7". This reduction applies also to neutral mesons such as 77, r]\ (f), cu according to the 
fraction of ss content in the meson itself. 

The hadron rates calculation proceeds through two steps; in the first one a primary number of 
hadrons emitted from the thermal source is calculated according to the formulae of previous 
section; in the second, all decays of these primary hadrons are performed according to known 
decay modes and branching ratios. The overall production of each particle is then computed by 
adding up the fraction stemming from decays of heavier particles to the primary one. The decay 
chain stops when tt, K, /j, or stable particles are reached, in order to match the numbers 
of production rates provided by experiments at e"^e~ colliders, including all decay products of 
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particles with cr < 10 cm. All hadrons included in JETSET 7.4 tables have been considered, 
with the quoted masses and widths. All other light flavoured resonances up to a mass of 1.7 
GeV have been included with masses, widths and branching ratios quoted by PDG [Q. The 
mass of resonances with F > 1 MeV has been distributed according to a relativistic Breit- 
Wigner function within ±2r from the central value. Experimental production rates used in the 



flt at ^/s = 91.2 GeV [10| and at = 29 ^ 35 GeV are weighted averages of all available 
measurements. 

The results of the flt are summarized in Tables 1, 2 and Figs. 2, 3. The errors on the parameters 
T, V and 7^ estimated from the flt include also uncertainties on input data such as branching 
ratios, masses and widths of particles involved in the flt. The flt procedure is described more 
in detail in ref. p. 

An excellent agreement between experimental and calculated production rate values is achieved. 
A signiflcative deviation is observed at y/s = 91.2 GeV for the S* production rate whose 
experimental value is still controversial because of the large disagreement between different 



experiments 11 



By using as input data the fltted T, V and 7^ values, the model also yields a deflnite prediction 
for the relative production rates of heavy flavoured hadrons once Rc = a{e~^e~ — > cc)/a{e~^e~ — *• 
hadrons) and Rb = a{e~^e~ bb)/a{e~^e~ — > hadrons) are known. 

In Table 3 the predicted multiplicities of some of the lightest heavy flavoured hadrons are 



quoted together with corresponding experimental vaues [|12|. A very good agreement with 
data is found. Two main sources of systematic uncertainties have been considered. The flrst 
one is simply a possible systematic error in the parameters flt which has been estimated by 
excluding the data points deviating the most from theoretical values and repeating the flt 
itself. The second one is related to possible different values of as a function of the primary 
quark mass. This has been taken into account by introducing two dimensionless variables 
Xc,Xb € [0, 1] such that XcV is the volume in e~^e~ cc and XbV in Z —>■ bb events, whereas 
V{1 — XcRc — XbRb)/{Ru + Rd + Rs) is the volume for Z ^ qq where g is a light quark. 
The measurement of the averaged charged multiplicity in e'^e^ —>■ cc and e~^e~ —>■ bb n 



determine x^ = 0.88 ± 0.03 and Xb = 0.70 ± 0.016 at ^ = 91.2 GeV and Xc = 0.89 ± 0.12 and 
Xb = 0.47 ± 0.15 at y/s = 29 -h 35 GeV. Inserting those values in the flt, it is found that new 
fltted parameters are almost unchanged. Other systematic effects related to uncertainties on 
Rg are found to be negligible. The variations of the fltted parameters due to those mentioned 
effects are summarized in Table 4. 

Another important issue is the independence of the fltted parameters, in particular T and 7^, 
from the neglected part of the hadronic mass spectrum, i.e. light-flavoured resonances with 
mass > 1.7 GeV. This has been tested by moving the cut-off point down to 1.3 GeV by steps of 
0.1 GeV. The new values are found to be consistent with the main quoted values within the flt 
error, as shown in Fig. 4. Morevoer, the number of primary hadrons stabilizes around 17 as the 
cut-off points increases. This indicates that missing heavy resonances production is negligible 
and should not affect signiflcantly the flt. 

One of the most relevant results is that temperatures values at both centre of mass energies 
turn out to be very close. If we add in quadrature the systematic uncertainties quoted in Table 
4 to the flts error quoted in Table 1 we get: 

Tlep = 162.8 ± 2.9 (12) 
TpEP = 169.5 ±3.7. (13) 
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The of the weighted average of those values is 1.8, which indicates a good degree of com- 
patibility. Also, it should be taken into account that the temperature value of PEP-PETRA 
might be affected by an additional systematic error due to having averaged production rates 
measurements spread between 29 and 35 GeV. 

4 Conclusions 

The problem of the relative production of hadrons in e^e~ collisions is treated with a thermo- 
dynamical approach postulating that hadronic jets must be identified with hadron gas phases in 
thermodynamical equilibrium before the hadrons themselves decouple and subsequently decay 
into lighter particles. Nevertheless, since the particle production depends almost linearly on 
the volume, the occurrence of additional jets, i.e. phases, in the same event should affect the 
relative hadron production rates very mildly provided that the temperature of the additional 
jets is the same. 

It has been shown that this model is able to fit very well the average multiplicities of light 
hadrons per hadronic event observed at LEP (^5=91-2 GeV) and PEP-PETRA (^s = 29^35 
GeV) for a variation of production rate by four orders of magnitude, from tt to fi. 
The predicted rates of heavy flavoured hadrons are also in good agreement with data. Only 
three parameters are required in order to reproduce data correctly, namely the temperature 
and the volume of the hadron gas and a parameter 7^ G [0, 1] allowing an incomplete strange 
chemical equilibrium. The temperature values at ^/s = 91.2 GeV and ^/s = 29 -j- 35 GeV, de- 
termined with the fit, are very similar and indicate a constant hadronization temperature. The 
thermalization of the system could be a characteristic of the quark-hadron transition, brought 
about by strong interactions. 
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Figure captions 



Figure 1 Behaviour of the jet partition function Z as a function of baryon number N for [S, C, B) = 
(left) and as a function of the strangeness S for (N, C,B) = (right), for T — 163 MeV 
and V — 21 Fm^. The suppression factors Z{Q — qj)/Z(Q) are smaller for increasing N 
and S. 

Figure 2 Fit of hadron production measured at LEP. Above: the solid line connects fitted values 
while data are shown as black dots. Below: fluctuations of measured points from the fitted 
values (solid line) in standard deviation units. Error bars include also contributions from 
uncertainty on masses, widths and branching ratios of hadrons. 

Figure 3 Fit of hadron production measured at PEP-PETRA. Above: the solid line connects fit- 
ted values while data are shown as black dots. Below: fluctuations of measured points 
from the fitted values (solid line) in standard deviation units. Error bars include also 
contributions from uncertainty on masses, widths and branching ratios of hadrons. 

Figure 4 Fitted temperature T, volume V, 7^ and number of primary hadrons at y/s — 91.2 GeV 
as a function of the light flavored resonance cut-off mass. 
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Tables 



Parameters 


v^ = 91.2 GeV 


v/i = 29 ^ 35 GeV 


Temperature(MeV) 


162.9±2.1 


169.3±3.5 


Volume (Fm^) 


21.4±1.9 


9.3±1.4 




0.696±0.027 


0.811±0.046 


xVdof 


38.6/17 


29.3/12 



Table 1: Fit results. 
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yi = 91.2 GeV 


v/i = 29 35 GeV 


Hadrons 


Calculated 


Measured 


Calculated 


Measured 


7r° 


10.03 


9.19 ±0.73 ±0.37 


6.186 


5.6 ±0.3 ±0.22 


7r+ 


8.952 


8.53 ±0.22 ±0.34 


5.415 


5.15 ±0.2 ±0.21 


K+ 


1.065 


1.18 ±0.064 ±0.013 


0.750 


0.74 ±0.045 ±0.008 




1.010 


1.006 ±0.016 ±0.013 


0.704 


0.74 ±0.035 ±0.008 


'n 


0.866 


0.94 ±0.11 ±0.11 


0.535 


0.61 ±0.07 ±0.0056 




0.939 


1.11 ±0.14 ±0.011 


0.544 




p' 


1.186 


1.29 ±0.13 ±0.17 


0.704 


0.81 ±0.08 ±0.11 


K*+ 


0.350 


0.357 ±0.034 ±0.089 


0.246 


0.32 ±0.025 ±0.006 




0.344 


0.379 ±0.020 ±0.089 


0.241 


0.28 ±0.03 ±0.006 


P 


0.534 


0.489 ±0.05 ±0.041 


0.293 


0.32 ±0.025 ±0.015 


v' 


0.0993 


0.22 ±0.07 ±0.0002 


0.0702 


0.26 ±0.1 ±0.00007 


</• 


0.123 


0.107 ±0.009 ±0.002 


0.106 


0.085 ±0.011 ±0.002 


A 


0.1(32 


0.185 ± 0.007 ± 0.007 


0.101 


0.103 ± 0.005 ±0.005 


(E+ + E-)/2 


0.0366 


0.044 ±0.007 ±0.0035 


0.0226 




SO 


0.0383 


0.036 ±0.007 ±0.0035 


0.0243 




A++ 


0.0925 


0.062 ±0.032 ±0.011 


0.0488 






0.0122 


0.0128 ±0.0007 ±0.0004 


0.00893 


0.0088 ±0.0014 ±0.0003 


(S+* + E*-)/2 


0.0192 


0.011 ±0.002 ±0.0014 


0.0111 


0.0085 ± 0.002 ± 0.0015 




0.00421 


0.0031 ± 0.0006 ± 0.00007 


0.00279 






0.000856 


0.00080 ±0.00025 ±0.0001 


0.000659 


0.007 ±0.0035 ±0.00015 



Table 2: Measured hadrons average multiplicities per hadronic event compared to the fitted 
values. The first error is the experimental one, the second error is due to uncertainty on hadron 
masses, widths and branching ratios. 
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= 91.2 GeV 


v/i = 29 H- 35 GeV 


Hadrons 


Calculated 


Measured 


Calculated 


Measured 




0.231 


0.221 ± 0.012 


0.256 


0.225 ± 0.035 




0.0914 


0.087 ±0.008 


0.106 


0.085 ±0.015 


D*+ 


0.107 


0.088 ± 0.0056 


0.113 


0.192 ±0.036 


B 


0.0902 


0.097 ±0.026 


0.0387 






0.694 


0.745 ± 0.069 


0.697 






0.0311 


0.037 ±0.009 


0.0339 


0.055 ±0.025 


Ds 


0.0536 


0.041 ± 0.0076 


0.0542 





Table 3: Measured heavy flavoured hadrons average multiplicities per hadronic event compared 
to predictions of the model. The experimental values have been averaged according to the 
procedure of ref. Values of i?^ = 0.17, Rb = 0.22 at = 91.2 and of Rc = 0.36, Rb = 0.09 
at ^/s = 29 -^ 35 have been used. 



Parameter 


v/i = 91.2 GeV 


v/i = 29 ^ 35 GeV 


Temperature(MeV) 


2.06 


1.08 


Volume(Fm^) 


1.34 


0.47 


Is 


0.017 


0.018 



Table 4: Systematic uncertainties on fltted parameters. 
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T = 163 MeV 



V = 21 Fm^ 
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